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Abstract 

We present a phase-space representation of quantum state vectors for two- 
cluster systems. Density distributions in the Fock-Bargmann space are con¬ 
structed for bound and resonance states of ®’^Li and ^’®Be, provided that all 
these nuclei are treated within a microscopic two-cluster model. The density 
distribution in the phase space is compared with those in the coordinate and 
momentum representations. Bound states realize themselves in a compact 
area of the phase space, as also do narrow resonance states. We establish 
the quantitative boundaries of this region in the phase space for the nuclei 
under consideration. Quantum trajectories are demonstrated to approach 
their classical limit with increasing energy. 
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1. Introduction 

The idea of formulation of quantum mechanics in a phase space is dis¬ 
cussed in numerous theoretical papers [D El El la El E]. The majority of such 
investigations are concentrated on establishing the link between quantum 
and classical mechanics. Due to the uncertainty principle there is no unique 
dehnition of the phase space. By this reason, different quantum phase space 
distribution have been proposed. In particular, in the phase-space represen¬ 
tation of Wigner and Husimi a quantum state is represented by a distribution 
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function (see the definitions, for instance, in PQ[2]), and the equations of mo¬ 
tion are of the Liouville type. 

A state-vector representation is another possibility to describe the dy¬ 
namics of a quantum system in phase space. In this case a quantum state 
is represented by a wave function, and the equations of motion are of the 
Schrodinger type. The dehnition of the phase-space representation is related 
to the choice of an operator which should be diagonal in this representation. 
In the coordinate representation the coordinate operator is diagonal, while 
the momentum operator is non-local. At the same time, the momentum 
representation diagonalizes the momentum operator, and the coordinate op¬ 
erator is non-local. Obviously, the coordinate operator and the momentum 
operator can not be diagonal simultaneously due to the uncertainty principle. 
Hence, one should seek for another operator. 

In Ref. [3] , mention was made that representation of a quantum state as 
a probability amplitude depending on two real variables related to the coor¬ 
dinate and momentum dates back to the papers of Fock [7] and Bargmann 
[8]. In the Fock-Bargmann space, a quantum state is represented as an entire 
function of a complex variable, with real and imaginary part of this variable 
being proportional to the coordinate and momentum, correspondingly. The 
Fock-Bargmann representation diagonalizes the creation operator, while the 
annihilation operator is non-local. 

A quantum problem can be resolved in any one of the above-mentioned 
representations. The Fock-Bargmann image of a wave function can be ob¬ 
tained from the wave function in the coordinate representation by a linear 
mapping, while the Husimi and the Wigner distribution functions are bilinear 
with respect to the wave function in the coordinate representation. However, 
the Fock-Bargmann representation is closely related to the Husimi distribu¬ 
tion. The latter distribution is equal to the square of the Fock-Bargmann 
image of the corresponding wave function multiplied by the Bargmann mea¬ 
sure. 

In Refs. [HE], Torres-Vega and Frederic suggested a quantum-state vec¬ 
tor phase-space representation. The authors postulated the existence of a 
complete basis of states \q,p) such that in the phase space a quantum state 
1-0) is represented by an £^(2) wave function '^(p, q) = (g,p|'0)- Here q,p are 
real values and the operators of coordinate and momentum in this basis take 
the form: 
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Thus Torres-Vega and Frederic developed a wave-function formulation of 
quantum mechanics in a phase space. Wave functions are governed by the 
Schrodinger equation in the phase space, while square of absolute value of 
the wave function plays the role of the probability density in the phase space 
obeying the Liouville equation. However, the quantum-state vector phase- 
space representation is not uniquely dehned, because there exists an inhnite 
number of bases depending on two real variables p, q which result in the 
foregoing expression of the coordinate and momentum operators Q and P. 

In Ref. [5] , the authors used coherent states as basis vectors and demon¬ 
strated that any coherent state representation leads to expression Q for the 
coordinate and momentum operators Q and P. Moreover, they concluded 
that only the coherent state representation makes it possible to dehne the 
operators of coordinate and momentum in such a manner. 

Following Klauder and Perelomov, in [5] a set of coherent states is defined 
as a result of action of the Weyl operator D{q,p) (a translation operator in 
the phase space) 

,p) =exp|^ {pQ-qP^ 
to any normalized vector |x) : 

\q,p-,x) = D{q,p)\x)- 

For any hxed vector |x) the set of coherent states \q, p; x) ensures a continuous 
representation of quantum states where the expansion coefficients can be 
interpreted as the wave function in the phase space. 

Then ^ ’’dilute” probability density in the phase space that 

equals the probability for a system to be localized in some ’’dilute” neigh¬ 
borhood of the center of the displaced state y. The degree of ’’diffusiveness” 
depends on the choice of vector y. If we choose the ground state of a har¬ 
monic oscillator as the vector x then |'0x(5') p) P is probability to hnd a system 
in the elementary phase volume AqAp = h near the point {q,p)- 

The Schrodinger equation in the phase space does not depend explic¬ 
itly on the vector x- However, in calculating mean values of operators the 
wave functions used for averaging should belong to the same coherent state 
representation, i.e., to the same vector x- This could be ensured by the re¬ 
quirement for the vector x to be an eigenvector of a certain operator. Hence, 
to formulate quantum mechanics unambiguously, one should solve two equa¬ 
tions in the phase space. 


D{q 
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It is possible to formulate quantum mechanics in a phase space so that 
one equation could be sufficient both for the solution of a quantum problem 
and for specihcation of the representation. As concluded in [5], the Fock- 
Bargmann representation suggests the best answer to this question. The 
Fock-Bargmann representation is a state-vector representation in the complex 
plane which is based on the mapping of the pair of bosonic creation and 
annihilation operators (a,a^) : 

a —)■ —, —)■ R. (2) 

In this representation quantum mechanics can be uniquely determined due 
to the fact that the creation operator is diagonal in this representation. 
The relation [a, a^] = 1 is followed by the relationships between the bosonic 
operators and the coordinate and momentum operators Q and P: 




a 




where b is the oscillator length. 

For a specihc hamiltonian H{Q,P) one can £x the relationship between 
the operators Q, P and bosonic operators (a, a'*') by choosing the value of the 
oscillator length b. Then the Hamiltonian can be expressed in terms of the 
bosonic operators H{Q, P) —)■ P(a,a^). Hence, the vector y is completely 
determined by the value of b. In this case the vector y is the vacuum vector 
| 0 >. 

The Fock-Bargmann representation can be associated with any Glauber 
coherent state representation based on the eigenvector y of the annihilation 
operator a. According to Perelomov [9], a coherent state describes a non- 
spreadable wave packet for an oscillator. Besides, coherent states minimize 
the Heisenberg uncertainty relation AqAp > h. Hence the coherent states 
are quantum states which resemble classic states the most. The Schrodinger 
equation in the Fock-Bargmann representation has the form: 


o o / ^ \ 

ih—IV’) = 'H(a,a^)|V’) -t ih—'tpFB{R) = P ( 'iPfb{R), (3) 
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Hence, a complete basis of the wave functions belonging to the class of co¬ 
herent states which minimize the uncertainty relation can be unambiguously 
determined as a result of solution of the Schrodinger equation in the Fock- 
Bargmann representation ([^. 

We follow precisely this strategy and present the wave functions and the 
probability distributions in the Fock-Bargmann representation. However, we 
don’t solve an equation of the type (|^, we use other way for obtaining an 
exact wave function of a quantum system. We concentrate our attention 
on the analysis of nuclear systems with a pronounced two-cluster structure, 
whereas authors of Refs. PE] have studied one-dimensional systems wherein 
the solutions could be found analytically. Much prominence is given to the 
determination of those regions of phase space, which are the most impor¬ 
tant for the dynamics of two-cluster systems. We also suggested the way 
of investigation of the density distributions in the phase space in a three- 
dimensional case, when the probability density distribution depends on six 
variables: absolute values of coordinate and momentum and four angles. 

There are some similarities between our approach and the Antisymmetrized 
Molecular Dynamics (AMD) [Tnl ITT] and the Fermionic Molecular Dynam¬ 
ics (FMD) [121 ESI- The AMD and FMD have been intensively employed to 
study a cluster structure of atomic nuclei and both of them appeal to a phase 
space. All three methods make use the same single particle orbitals to con¬ 
struct a many-particle wave function in the form of the Slater determinant. 
It means that all methods involve the same part of the total Hilbert space, 
provided that all of them take into account the same partition (or cluster¬ 
ization) of A nucleon system. However, there are more differences between 
our approach and AMD and FMD. First, we use the Slater determinants 
as the generation functions for a complete basis of many-particle oscillator 
functions, describing the most important, from physical point of view, types 
of motion of many-particle system. In our approach cluster parameters are 
the generator coordinates that allow us to select necessary basis functions 
from the inhnite set of oscillator functions. Meanwhile, in the AMD and 
FMD they are independent variables in the phase space. By using a set of 
many-particle oscillator functions, we reduce the Schrodinger equation to the 
matrix form. The AMD and FMD make use of the time-dependent equations 
derived from the time-dependent variational principle. Having obtained the 
wave function in a discrete, oscillator representation, we then transform it to 
the Fock-Bargmann or phase space, where we analyze phase trajectories of 
a quantum many-particle system. 
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It should be emphasized that all the results presented in this paper are 
obtained within the phase-space formulation of quantum mechanics which 
is valid both for hnite h and when h goes to zero. We consider some sim¬ 
ple model cases as well as two-cluster systems. The model cases, such as 
a three-dimensional harmonic oscillator and a plane wave, help us to reveal 
peculiarities of phase space portraits. However, our main aim is to study 
phase portraits of real physical systems, namely, light atomic nuclei. All cal¬ 
culations are performed within a microscopic two-cluster model based on the 
resonating-group method (RGM) [T3j. Our analysis started from assump¬ 
tions of the following cluster structure of the nuclei under consideration: 

= a + d, Li = a+ ^Be = a+^ He, ^Be = a + a. Instead of solving an 

integro-differential equation in a phase space as the authors of Ref. |lj do, we 
deal with the RGM Hamiltonian in the representation of the Pauli-allowed 
harmonic-oscillator states dehned in the Fock-Bargmann space. Doing so, 
we reduce the integro-differential equation to a set of linear equations for 
coefficients of the expansion of the wave function in the harmonic-oscillator 
basis. 

Our paper is organized as follows. In Section]^ we introduce all necessary 
dehnitions and formulate an approach for constructing density distributions 
in the Fock-Bargmann space. In Section the effectiveness of the suggested 
approach is demonstrated for two simple model problems: a harmonic oscil¬ 
lator and a free motion in 3D space. In Section]^ we give a brief review of the 
employed two-cluster model. Details of the calculations are shown in Section 
1^ Phase portraits for bound and resonance states in the light atomic nuclei 
are presented in Section]^ Finally, in Sectionwe conclude the main results 
obtained. 

2. The Fock-Bargmann representation 

In this section and Section we shall use dimensionless units for energy 
and length. The energy E is measured in units of and the length is 
measured in terms of oscillator length b. The value of the oscillator length 
will be determined in what follows. 

The transition from the wave function in the coordinate or momentum 
representation to the wave function in the Fock-Bargmann representation is 
performed using the Bargmann-Segal integral transformation 



( 4 ) 
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the kernel of which is the modihed Bloch-Brink orbitals |15j : 

r) = ^ exp + V2 (rR) - . 

In the Fock-Bargmann space, wave functions are entire analytical functions 
of complex variable R, while ^ and ff are coordinate and momentum vectors, 
respectively. 

The modihed Bloch—Brink orbital is an eigenfunction of the coordinate 
operator in the Fock-Bargmann space: 

r) = riF(R,r), ? = ^ (R + 

On the other hand, the modihed Bloch-Brink orbital is a coherent state and 
generates a complete basis of functions of a 3D-harmonic oscillator. 

Having calculated the wave function in the Fock-Bargmann space, we can 
hnd the density distribution, which depends on six variables: absolute values 
of the coordinate and momentum rj and four angles Qff: 






Ti5(R)|2exp 



d^dff 
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( 6 ) 


Here 
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exp{ — (R ■ S)} 


d^dff 
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R* 


is the Bargmann measure. 

In essence, dDE{^, if) is a density matrix in the Fock-Bargmann space. It 
is positively dehned for all values of ^ and ff as opposed to the Wigner func¬ 
tion, which is often used to construct a density distribution in a phase space. 
In fact, rf) coincides with the Husimi distribution within dehnition of 

variables ^ and ff. 

It is rather difficult to investigate the density distribution DE{^,ff), be¬ 
cause it depends on six variables. To solve this problem and extract as much 
as possible physical information, we propose to integrate DE{^,ff) over the 
solid angles and Qff. Then we come out with the density distribution 
De{^,ti) which depends only on the absolute values of the coordinate and 
momentum. 


De{^,v) 


f dfl^ f dQffdDE{i,ff) 
d^di] 
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It is important to recall, that the variable ^ is analogue of the distance 
between interacting particles (clusters) and the variable p represents a mo¬ 
mentum of relative motion of particles (clusters). 

Given the energy E and other integrals of motion, like the orbital mo¬ 
mentum and parity, the density distribution De{^,ti) is comprised of the 
inhnite number of phase trajectories. The phase trajectories are determined 
as a continuous set of points in the (^, rj) plane for the hxed values of the 
density distribution = const, while the probability of realization of 

the phase trajectory is proportional to the values of De{^,ti). Hereinafter 
we shall call an inhnite set of phase trajectories characterizing a particular 
quantum state of the system under study the phase portrait. 

With increase in energy E, all the quantum phase trajectories gradually 
approach their classical limit. Therefore, analyzing phase portraits of quan¬ 
tum systems we can make a quantitative estimate of the energy such that the 
maximum of the density distribution De{^, f]) falls on the classical trajectory. 

In [15], we have constructed and analyzed the phase portraits for the free 
motion of a ID quantum particle and for the motion in the held of Gaussian 
potential. In the present paper we concentrate on the three-dimensional case. 
Before proceeding to the analysis of phase portraits of two-cluster systems, 
let us discuss the density distribution for the states of a three-dimensional 
oscillator with the number of quanta N and for a free motion of a 3D-particle 
in the states with the orbital momentum 1. 

3. Model problems 

3.1. Harmonic oscillator 

The overlap integral of the modihed Bloch-Brink orbitals 

(S|R) = j i^(S, r)it'(R, r)(ir = exp(RS) 

generates a complete basis for the harmonic-oscillator functions in the Fock- 
Bargmann representation: 

OO 

exp(RS) = 5^ 5^ ^ 

n=0 I 
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Here both the absolute value and the solid angle of the complex vector R 
take complex values: 


R=V^ 



p2)+2^(fp), 


^R=- 

" R 


-V^) + 2*(fh) 


Each basis function corresponds to the total number of quanta N = 2n+l, 
where I is the orbital momentum and n is the number of radial quanta: 




An 


{2n)\\i2n + 21 + l)\\ 


( 8 ) 


Hence, according to the dehnition ([^, a density distribution for the |n, /; R) 
state of harmonic oscillator should be written as 


dDn,i{^,rf) = 


21 + 1 


|n, /, m; S)(n, /, rri] R| d^B- 


Here averaging over projection m of the orbital momentum is performed. 

Integrating over the solid angles related to the coordinate ^ and momen¬ 
tum f], we come to the expression 


DnAU) = IS = 


(2n) 
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where 


1 3 

2 ’ 2 ’ ie-vTJ' 

I {-in2l-2u)\ 
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As observed in the right panel of Fig. with increase in the orbital momen¬ 
tum I the density distribution Dn,i{^, v) is slightly forced out from the region 
of small ^ and rj and narrows. 

The left panel of Fig. [T] demonstrates that the density distributions 
Dn^i^, v) corresponding to the same total number of quanta N, but different 
orbital momenta I, peak on the same circle: 

e+v^ 


= 2n + I + 2. 


(9) 
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Eq. ([^ describes a classical trajectory of the particle with the energy = 








Figure 1: Phase portraits of a 3D oscillator = const with the total num¬ 

ber of quanta 2n -I- ? = 20 as a function of dimensionless coordinate ^ and momentum 
rj. The left panel corresponds to Dn=io,i=o{^,r]) = const, the right panel demonstrates 
Dn= 5 ,i=io{^,r]) = const. The heavy line on both panels denotes the classical trajectory 
defined by Eq. (§. 

2?7, + Z + 2. Hence, the most probable locus of points in the phase space 
for the eigenstate of a 3D harmonic oscillator with the number of quanta 
N = 2n + I corresponds to the classical energy = N + 2 rather than to 
Ed = N + 3/2. The maximum of density distribution for the iVth state of a 
ID harmonic oscillator was matched by the classical energy E^ = N instead 
of Ed = N + 1/2 (see [U ITH]). 

3.2. Plane wave 

In the Fock-Bargmann representation, a plane wave corresponding to the 
momentum k becomes: 

0k(R) = 7r"^/'^exp 1“ Y “ iA/2(kR) + ^| • (10) 

It is an eigenfunction of the momentum operator k: 

k0k(R) = k0k(R), k = -^ (R - Vr) . 
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The density distribution for the plane wave (10) has the Gaussian dependence 
on the momentum rj and is unaffected by the coordinate ^ : 


= |0k(R)pd/is = TT ^/^exp (-(k-p)2) (11) 

Obviously, the density distribution ( pT| peaks on the line rf = k. This line 
coincides with the classical phase trajectory of a free particle with momentum 

k. 


The problem becomes more intricate if we consider free motion of a par¬ 
ticle in the states with given value of the orbital momentum 1. We can 
expand the eigenfunction (10) of the momentum operator k in the basis of 
the harmonic oscillator functions ([^ in the Fock-Bargmann space: 


|n,(12) 

n=0 / m 


Expansion (12) in the Fock-Bargmann representation is much the same as 
the multipole expansion of the plane wave in the coordinate representation. 

The expansion coefficients Cni{k) coincide with the basis functions of 
harmonic oscillator in the momentum representation. In particular, for zero 
orbital momentum the expansion coefficients Cni{k) are defined as 


(-y) ■ (13) 

Generally in the Fock-Bargmann space any wave function T^(R) of a quan¬ 
tum state characterized by the orbital momentum I, its projection m and 
energy E can be represented as the expansion into the basis of the harmonic 
oscillator functions ([^: 


n=0 I 


where the expansion coefficients are the solutions of a set of linear alge¬ 
braic equations to which the Schrodinger equation is reduced. The expansion 
(14) is written in the Fock-Bargmann space; it can be also presented in coor¬ 
dinate and momentum spaces with the same set of the expansion coefficients 
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{C^}. Such way of representing and calculating the wave functions is the 
essence of the method we follow in the present paper. 

Hence, the density distribution for the quantum state described by the 
wave function takes the form: 


dD^Al v) = ^ E E E C5Ci,\k, I, m; S)(k, I, m; R| dfig. 

n=0 n=0 m 


Upon integrating the foregoing expression over the solid angles and 
we obtain: 




e+v^ 
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n=0 n=0 
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E< 


u=0 

X Im 


e+v^ 


2 ^2\n+n+2zx 


(r - v^) 


, 1 fri 

H ( - + ^ ^2 _ ^ 2 ’^ + u + l,n + u + 1 


where B{z-, a, b) is the incomplete beta function. 

Substituting the coefficients ( |I^ in the formula given above we derive 
the density distribution for a free motion of the particle with zero orbital 
momentum. Figure [^demonstrates the phase portraits for a free particle with 
k = 1 and /c = 5. As observed in Fig. [^ for k = 1 the density distribution 
oscillates in the variable These oscillations, resulting in a number of closed 
trajectories, are of quantum nature. With increasing the energy the density 
distribution becomes smoother; all the quantum trajectories become inhnite 
and approach the classical trajectory r) za k. 

Figure [^compares the density distributions in the Fock-Bargmann space, 
coordinate representation and momentum representation for a free 3D par¬ 
ticle with the energy E = fc^/2, k = 1 and zero orbital momentum. 

In the coordinate representation the wave function of a free particle in 
the state with / = 0 is described by the spherical Bessel function ipk (r) = 
^^/2/7rjo{kr), while in the momentum representation it is just the Dirac delta- 
function 'ijjk (p) = 5{p — k). As is clear from Fig. [^ locations of the maxima 
for the density distributions in the variable ^ in the Fock-Bargmann and 
coordinate representations almost coincide, while the density distribution in 
the variable r] peaks at r; ~ 1.25 instead of r/ = 1 as the momentum density 
distribution does. 
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Figure 2: Phase portraits of a free 3D particle with the energy E = fc^/2 and zero orbital 
momentum as a function of the dimensionless coordinate ^ and momentum rj. The left 
panel corresponds to A: = 1, while the right panel corresponds to A: = 5. 


4. Two-cluster model of light nuclei 

We present shortly main ideas of the two-cluster model that we are going 
to use to study the dynamics of two-cluster systems. First of all, we restrict 
ourselves with the lightest nuclei of p-shell. 

A wave function for the two-cluster partition A = Ai + A 2 is represented 
as 

= d{|4>. (Ai) $2 (A,)],i.is («) n (q)} j. (15) 

where A is the antisymmetrization operator, q is the Jacobi vector which is 
proportional to the vector r connecting the centers of mass of the interacting 
clusters 


q = r 


A 1 A 2 


Ai -\- A2 V Ai J- A2 


A1A2 


Ai " 




iezli 


jeA2 


(16) 


It is assumed that the number of nucleons in each cluster Ai does not exceed 
4: 1< Ai <4 . As one sees, we use the LS scheme of coupling, when the 
total spin S is coupled with the total orbital momentum L and generates the 
total angular momentum J. For two interacting s-clusters, the total orbital 
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Figure 3: The left panel presents the density distributions DE,i=o{imax,v) the phase 
space (a solid line) and De,i=o{p) in the momentum space (a dashed line) for a free 
3D particle with the energy E = k'^/2, k = 1 and zero orbital momentum. The right 
panel shows the density distributions Dej=o{^, Vmax) in the phase space (a solid line) and 
De,i=o(p) in the coordinate space (a dashed line). 


momentum L coincides with the angular momentum of the relative motion 
of clusters. 

To hnd the wave function (q) of relative motion of the clusters one 
has to solve a set of integro-differential equations (see the details in Ref. 
HZ!)- The set of equations and the form of the wave function ( [l^ are the 
key elements of the well known Resonating Group Method. 

An integral part of the integro-differential equations for the functions 
'4’ls (^) antisymmetrization operator A and originates from 

the potential, and kinetic energy operators, and the norm kernel. The latter 
means that interaction between clusters is energy dependent. By neglecting 
the operator A, we obtain differential equations with a local potential. This 
approximation is called the folding approximation or the folding model. 

The Form and procedure of solving the equations for the functions (q) 
can be simplihed, if we use a set of square-integrable functions. Such way 
of solving the equations of the Resonating Group Method with oscillator 
functions is called an Algebraic Version of the Resonating Group Method [T8l 
To realize the algebraic version, we introduce a set of cluster oscillator 
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functions 

\n, L, S, J) =!{[$! (Ai) $2 (/l 2)]5 <^>nL (?) (q)}^ , (17) 

where 


*hnL (r, h) 
P 


{-iTMnL (/), 

r _ / 2r (n + 1) 

6’ VT{n + L + 3/2) 


(18) 


is an oscillator function in the coordinate space {b is the oscillator length). 
We are interested in wave functions and density distributions both in the 
coordinate and in momentum space. They will be constructed with the help 
of the oscillator functions in the momentum space 

P = P-b, 


The cluster oscillator functions \n, L, S, J) are totally antisymmetric and 
compose a complete set of basis functions with specific physical properties. 
Namely, they are a part of the total Hilbert space describing the Ai + A 2 
clusterization of the system of A nucleons with fixed internal clusters func¬ 
tions $1 (Hi) and $ 2 (^ 2 ). The cluster oscillator functions can be used to 


expand any wave function of the type (15): 




( 20 ) 


LS n=0 


The expansion coefficients {Cni} represent a two-cluster wave function in a 
discrete, oscillator representation and obey the system of linear equations 


EE 

LS m=0 


n, L, S, J 


H 


m,L,S, J) — E ■ {n, L, S, J\m, L, S, J) 


CmL = 0 , 


where 


n, L, S, J 


H 


m, L, S, J 


( 21 ) 

is a matrix element of the hamiltonian be¬ 


tween the cluster oscillator functions and (n, L, S', J|m, L, S', J) is an overlap 
of these functions or the norm kernel. For two s-clusters, the norm kernel 
has a very simple form 


(u, Ly S', Ly S', (/) bfi,m^n- 


( 22 ) 
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The constants A„ are the eigenvalnes of the antisymmetrization operator. 
The states In, L, S, J) with A„ = 0 are called the Panli-forbidden states. The 
antisymmetric basis fnnctions with the qnantnm nnmbers n, L, 5, J do not 
participate in constrncting the wave fnnction (15). Only the Panli-allowed 
states (i.e., the states with A„ > 0) take part in describing the dynamics of 
the two-clnster system. 


Eqnation (22) indicates that the clnster oscillator fnnctions 0 are not 
normalized to nnity despite the fact that the fnnctions <hi (^i), $2 (^ 2 ), 
^nL {q) and Ylm (q) are properly normalized. The antisymmetrization oper¬ 
ator A is responsible for that. By renormalizing the basis fnnctions 


\n 


,L,S, J) = \n,L,S, J)/VK 


and the expansion coefficients CnL = we arrive at the standard 

matrix form of the Schrodinger eqnation with the orthonormal basis 


LS m 


H 


?7Z, Z/, J J E ' Sn^m CfyiL 0, 


(23) 


where n and m nnmerates only the Panli-allowed states. 


Formally, the expansion (20) of the wave fnnction contains an infinite set 


of basis fnnctions; however, actnally we need a large bnt restricted set of fnnc¬ 
tions. In an oscillator representation, sitnation is similar to the coordinate 
form of the Schrodinger eqnation where one needs to find a wave fnnction np 
to a finite distance Ra- Beyond this point the well-known asymptotic form 
of the wave fnnction is valid. Ra determine a distance where a short-range 
interaction is negligibly small and an asymptotic part of the hamiltonian is 
dominant. The same is trne for the discrete representation. We need to cal- 
cnlate the wave fnnction np to the finite valne of n = A'A; starting from this 
nnmber of qnanta an asymptotic form for the expansion coefficients of the 
wave fnnction is valid. Like Ra-, the parameter Na draw a border between 
an internal and asymptotic regions. When solving the Schrodinger eqnation 
nnmerically both in the coordinate and oscillator representations, the param¬ 
eters Ra and Na are used as variational parameters. One has to determine 
minimum values of Ra and Na such that further increasing of them do not 
change the results of calculations. 


To solve the system of equations (21) or (23), one needs to take into 


account the boundary conditions. The asymptotic form for a wave function 
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of a bound state in the coordinate space 






exp {-Kg}/g, 


2m \E\ 


K = 




as was shown in Ref. [19], is transformed into 


Rn — \/47i + 2L + 3 


(24) 


(25) 


for the expansion coefficients of the wave function. Similar relations are 
valid for the wave function of a continuous spectrum state {single channel 
case) in the coordinate space 

^LS (?) ~ sin (^kq + 

2mE 

n? 

and in the oscillator representation 

CnL ~ sin +(27) 




where hi is a phase shift. To avoid bulky formulae and long additional 


explanations, in equations (24)-(27) we showed an asymptotic form of the 
wave functions for neutral clusters. Similar expression can be written for 
charged clusters. 

There is an equivalent form for the expansion (20). We can write expan¬ 


sion for the intercluster wave function by using the same set of the expansion 
coefficients 


i'is i‘l) = Yl CnL<*nL (f, 5) 


(28) 


n=0 


Similar form can be used to determine intercluster function (p) in mo¬ 
mentum space. The functions 'ijj'lg (g) and "0^5 {p) cire connected by the 
Fourier-Bessel transformation 


^LS iP) = dqq^jL (pq) (?) 


(29) 
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One important note should be made. The wave function (15) for bound 
states is traditionally normalized to unity 




= 1 


n=0 


but it is not the case for the corresponding intercluster functions 

{i’isli’is) = ^Lj- (30) 


In the oscillator representation it reads as 

OO OO 

SlJ = \CnLf/Xn. (31) 

n=0 n=0 

The deviation of the quantity Slj from unity shows how strong is the effect 
of the Pauli principle. 

Having calculated the expansion coefficients we can easily con¬ 

struct the wave functions and density distributions in the coordinate, mo¬ 
mentum, and the Fock-Bargmann or phase spaces. 

We do not dwell on the calculation of matrix elements of the hamiltonian 
between the cluster oscillator functions (HIl' We refer the reader to the 
review [20], where all necessary formulae are presented. They could help 
one to calculate the matrix elements of the hamiltonian, which include the 
central and spin-orbital components of the nucleon-nucleon forces and the 
Coulomb interaction as well. 


5. Details of calculations 

Our main objective is to study light nuclei with a pronounced two-cluster 
structure. Among these nuclei are the '^Li and ^Be, because the 

two-cluster decay thresholds a + d, a II and a -1-^ He, respectively, he 
not far from the ground state of the nuclei and other two- and three-cluster 
thresholds are at higher energy. There are the strong grounds to believe that 
the channels a + d, a+^ H and He are responsible to a great extent for 
the structure of bound and low-lying resonance states in Li and Be 

nuclei. We also consider the ^Be as a -|- a conhguration, which generates a 
set of the rotational 0^, 2'^ and 4'^ resonance states. 
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In our calculations we make use of the Minnesota nucleon-nucleon poten¬ 
tial suggested by Tang and coworkers. The central part of the potential is 
taken from Ref. 1211 and the spin-orbital components are taken from Ref. 
[22] (IV version). 

In such type of calculations we have got two free parameters. The hrst 
parameter, the oscillator length b, we use to minimize the energy of the se¬ 
lected two-body threshold. In other words, we use this parameter to optimize 
description of the internal cluster structure. The second parameter u, deter¬ 
mining the odd components of the Minnesota potential, is htted to reproduce 
the bound state energies for all nuclei but ^Be. The latter nucleus has no 
bound state and the ’’ground state” actually is a very narrow resonance state 
(its energy i?=0.0918 MeV and width r=5.57 eV) that can be treated as a 
quasi-stationary state. Thus for ®i?e, we hnd the parameter u, which repro¬ 
duces fairly well the energy and width of the O’*" resonance state. To study 
peculiarities of the O’*" resonance state in the ^Be, we make two different cal¬ 
culations. One was mentioned above. This result we mark as a Resonance 
State (RS).In the second variant of calculation, which will be marked as a 
Bound State (BS), we switch off the Coulomb interaction and thus obtain a 
bound state in ®Re. 

To study effects of the Coulomb interaction in the mirror nuclei "^Li 
and ^Be, we use the same input parameters, which were adjusted for the 
ground state. In this case the difference in position of the bound and 
resonance states in '^Li and ^Be is due to the Coulomb interaction. 

In the present model, the total spin S and the total orbital momentum 
L are the quantum numbers. For the two-cluster conhgurations, which are 
taken into considerations, the total spin is determined by the second cluster, 
as the hrst cluster, alpha-particle, has the spin equals zero. 

We use the following scheme of calculations. First, we construct matrix 
elements of the hamiltonian and other operators of physical importance be¬ 
tween the cluster oscillator functions. We use Na= 200 oscillator functions in 
all our calculations. This number of basis functions provides us with conver¬ 
gent and stable results both for the bound states and scattering states as well. 
Second, we calculate the eigenvalues and eigenfunctions of the hamiltonian. 
As a result we obtain a bound state (if any) and a large set of pseudo-bound 
states. The latter are the states of continuous spectrum with wave functions 
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normalized to unity in a fixed basis of functions 


Na-l 

Y, |C„ (£„)!" = 1 

n=0 

and obeying the conditions 


CNa {Ea) — 0 , 

where is the energy of the ath (a =0, 1,2,..., Na-l) pseudo-bound state. 
By using the eigenfunctions we construct the density distributions 

in the coordinate, momentum, and Fock-Bargmann spaces. 

Third, we calculate the phase shifts of elastic scattering by solving a 
system of linear equations, which takes into account the proper boundary 
conditions for scattering states. It allows us to determine energy and width 
of resonance states. 

In Table we show the input parameters of the calculations and the spec¬ 
trum of bound and resonant states of the light atomic nuclei. We also indicate 
in the Table the dominant two-cluster channel taken into consideration. The 
experimental data are from Refs. [231121] • As we see, with the input param¬ 
eters, indicated in Table, we obtain a fairly good description of the bound 
and resonance states comparing to the experimental data. However, the en¬ 
ergy and width of few resonances slightly differ from the experimental values. 
This can be attributed to the peculiarities of the used nucleon-nucleon po¬ 
tential and restrictions of the present model (using, for instance, the same 
oscillator length for both interacting clusters). We consider this drawback of 
the present calculations not crucial for the interpretation and validity of the 
results which will be discussed bellow. 

6. Results and discussion 

We start our discussions with wave functions in the coordinate space. In 
Figure]^ we present the wave functions of the for 10 lowest O’*" eigenstates 

of two-cluster hamiltonian. One notices that all the wave functions have 
nodes approximately at the same points in the coordinate space. This feature 
is typical for all nuclei under consideration. Position and the number of nodes 
depends on the nucleus and its angular momentum. These hxed nodes results 
from the orthogonality of the wave function of intercluster motion to the 
wave functions of the forbidden states; and the number of nodes is equal to 
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Table 1: Spectrum of bound and resonance states of the light nuclei and the input param¬ 
eters of calculations. Calculated energy and width are given in MeV, the experimental 
values of energy and width are given in MeVikeV. The theoretical and experimental width 
of the O’*' resonance state in ^Be are presented in eV. 


System 

Input 

Theory 

Experiment 

Nucleus 

+ A 2 

b, fm 

u 


E 

r 

E 

r 

^Li 

a + d 

1.3110 

0.9254 

1+ 

-1.4750 

- 

-1.4743 

- 





3+ 

0.8480 

0.0284 

0.712 ± 2 

0.024 ± 2 





2+ 

4.2880 

3.0052 

2.838 ± 22 

1.30 ± 100 

’^Li 

a+^H 

1.3451 

0.969 

3- 

2 

-2.4676 

- 

-2.4670 

- 





1 - 

2 

-1.6040 

- 

-1.9894 

- 





7- 

2 

2.4710 

0.1285 

2.185 

0.069 





5- 

2 

4.9390 

1.7121 

4.137 

0.918 

'^Be 

a He 

1.3451 

0.969 

3- 

2 

-1.6302 

- 

-1.5866 

- 





1 - 

2 

-0.8161 

- 

-1.1575 

- 





7- 

2 

3.3360 

0.2232 

2.98 ± 50 

0.175 ± 7 





5- 

2 

5.7420 

2.0207 

5.14 ± 100 

1.2 

^Be 

cx cx 

1.3736 

0.950 

0+ 

0.0818 

2.40 

0.0918 

5.57 ± 0.25 





2+ 

1.2840 

0.6418 

3.12 ± 10 

1.513 ± 15 





4+ 

9.7970 

3.5827 

11.44 ± 150 

^ 3.500 


21 




the number of the Pauli-forbidden states. This fact is the key element of 
the Orthogonality Condition Model [251126], which is a simplihed version of 
the Resonating Group Method with an approximate treatment of the Pauli 
principle and a local cluster-cluster interaction. 

The wave functions in the momentum space have also nodes, but their 
structure is not so simple and pictorial. 




Figure 4: Coordinate wave functions of 10 lowest 0+ states in ®i?e. 

In Figure we compare the coordinate wave functions of the ground 
states in ^Li, "^Li, ^Be and ^Be (BS). The wave functions of and ^Be 
are indistinguishable in this Figure, despite that the energy difference of the 
bound states is 0.84 MeV. 

To demonstrate whether our basis of cluster functions is large enough to 
provide correct results, in Figure]^ we also display an asymptotic behavior 
of the wave functions. One can see that the wave functions are decreased as 
exp {—Kg} /g at large values of g. The order of the curves depends on the 
bound state energy: the smaller is the energy, the lower is the corresponding 
curve. 
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r, fm r, fm 


Figure 5: Wave functions of the ground state of '^Be and ®i?e in coordinate space. 

A short range behavior is shown on the left-hand side, while an asymptotic behavior is 
demonstrated on the right-hand side of the Figure. 
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6.1. Phase portrait of bound states 

Now we turn our attention to the phase portraits of bound states of the 
two-cluster systems. In Figures we display the phase portraits of the 

and ^Be bound states. The phase portrait for the ^Li ground state was 
shown in our previous publication [16]. The general feature of these figures is 
that the two-cluster system is concentrated in a rather narrow region of the 
phase space, despite the fact that some bound states are weakly bound ones. 
As was expected, the more dispersed is a state in the coordinate space, the 
more compact it is in the momentum space. And vise versa. It is interesting 
to note that a maximum of the density distributions in the phase space of 
the bound states lies at ~ 1.0 and 3.5 < < 4. The above-mentioned value 

of r] is noticeable differs from the dimensionless momentum k (which for the 
deepest bound state in Li equals fc=0.46). In our opinion, this indicates that 
the shape of the density distribution for a bound state has a pure quantum 
character. 

As we see from Fig. the wave functions and thus the density distribu¬ 
tions of the ground states in the coordinate space have nodes. There are also 
nodes in the momentum space for a wave function of bound states. However, 
the density distribution of the states in the phase space does not have any 
node in the range 0 < ?7 < cxd and 0 < ^ < cxd. 
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Figure 7: Phase portraits of the 3/2 (left) and 1/2 (right) bound states in ^Be. 

6.2. Phase portrait for resonance states. 

There are narrow and broad resonance states in the nuclei of interest. 
The narrowest resonance state is observed in ^Be. Indeed, the calculated 
width of the O’*" resonance state is only 2.40 eV. The broadest resonance 
state is also observed in ®i?e. The width of the d"*" resonance states exceeds 
3.5 MeV. One would expect quite different density distributions and shape 
of phase portraits for narrow and broad resonance states. This is so, as can 
be seen from Figures where we display the phase portraits for the 

resonance states in '^Li and ^Be. Each nucleus is represented by two 
resonance states, one of which is narrow (S’*" in ^Li, 11‘2~ in Li and '^Be) 
and the other is broad {2~^ in ^Li, 5/2“ in '^Li and "^Be). When we are saying 
”a narrow resonance state” or ”a broad resonance state”, we mean not only 
the absolute value of the resonance width, but also the value oiT/E. For 
narrow resonance states this ratio is below 0.067, and for broad resonance 
states it exceeds 0.34. For the narrowest O’*" resonance state in ^Be, the ratio 
T/E ^ 3 ■ 10“^ 

There is a strong resemblance between narrow resonance states and bound 
states. Both of them realize themselves in a compact area of the phase space. 

As for the broad resonance states, they have a principal maximum at 
relatively small values of coordinates f and momenta r]. Besides, they also 
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Figure 8: Phase portraits of the S'*" (left) and 2+ (right) resonance states in ^Li. 


have many regular maxima at a fixed value of t], but different values of 
This structure reflects an oscillatory behavior of the coordinate wave function 
and 5-like behavior of the wave function in the momentum space. 

In Figurejl^we compare phase portraits of the O’*" ’’ground” state in ^Be, 
obtained with and without the Coulomb interaction (RS and BS calcula¬ 
tions). The energy of the resonance state is 0.0818 MeV above the a + a 
threshold, while the energy of the bound state is -1.3529 MeV bellow the 
threshold. 

One can see, that the phase portraits of resonance and bound states 
are quite similar. This similarity can be attributed to the fact that in the 
internal region, where interaction between clusters is strong, wave functions 
have very close behavior. The asymptotic part of these functions are totally 


different. For a bound state, the function ^jJ (q) exponentially decreases (24), 


while it slowly decreases (26) for a resonance state. Thus, the asymptotic 


part of wave functions gives a small contribution to a density distribution 
in the Fock-Bargmann space. To prove that in the internal region wave 
function of resonance and bound state are similar, we display in Figures [T^ 
the density distributions in the coordinate and momentum representations. 
Indeed, these figures show that there is small difference between the wave 
functions or density distributions of the bound and resonance states in ^Be. 
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Figure 9: Phase portraits of the | 


and I resonance states in Li. 





Figure 10: Phase portraits of the | and | resonance states in Be. 
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Figure 11: Phase portraits of 0+ resonance (left) and bound (right) states in ®_Be. 


Maximum of density distributions is approximately in the same point of the 
coordinate space. However, maxima of density distribution of a BS and RS in 
the momentum space lie at different values of momentum p, which connected 
with different energies of the resonance state (0.0809 MeV) and bound state 
(-1.3529 MeV). Since the wave function of the bound state is more compact 
in the coordinate space than the resonance one, it is more dispersed in the 
momentum space. 

It is interesting to compare a density distribution in the phase space with 
density distributions in the coordinate and momentum spaces. For this 
aim, we hnd a point in the phase space (^max, hmax) such that the density 
distributions has maximum. And then, we calculate the density distributions 
-D ('Cl hmax) and H (Cmax,Thus from a large number of trajectories, we 
selected only two of them, the hrst one is with a hxed momentum pmax, 
the second one is with a hxed coordinate C- The hrst trajectory we compare 
with the density distribution in the coordinate space D (r), and the second 
trajectory is compared with the density distribution in the momentum space 
D (p). These results are presented in Figure [I^ One can see, that the density 
distributions D (C, Pmax) and D (r) are quite close to each other, they have 
maximum approximately at the same point ^ — r. There are more diherences 
in behavior of the density distributions D (Cmax,h) and D (p). Maximum of 
D (p) is shifted to smaller values of p, as compared to the function D (Cmax, v)- 
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r, fm p, fm 


Figure 12: Density distributions of the 0+ resonance and bound states in ®i?e in coordinate 
(left) and momentum (right) spaces. 
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Such relation between and D{r), -D(^max,h) D (p) is also 

observed for bound and resonance states of other nuclei. 




Figure 13: Density distributions in phase space and in coordinate and momentum repre¬ 
sentations for O"*" resonance state in ®Be. See more detail in text. 

By closing this section, in Figure [T4| we demonstrate the phase portraits 
for the broad 2~^ and 4+ resonance states. One can see that the oscillations of 
the density distribution along the axis ^ is more frequent for the 4'^ resonance 
state than for the 2'^ resonance state. This is expected feature of the density 
distribution, because the energy of the 4+ resonance state exceeds the energy 
of the 2+ resonance state by a factor of 8. Both resonance states reveal strong 
quantum effects, because they have prominent maximum in the region of 
small intercluster distances in the phase space. 

As with the bound states, a maximum of the density distribution in phase 
space for the resonance states is also observed at the values of variable p, 
which signihcantly differ from the corresponding value of the dimensionless 
parameter k. To demonstrate this, we consider the resonance states in ^Be. 
The dimensionless parameter k for ®i?e, where the 0+ resonance state has the 
smallest energy and the 4+ resonance state has the largest energy, varies at 
the range 0.086< k <0.944. However, a maximum of the density distribution 
for the 4+ resonance state is achieved at p ^ 1.5. This value has to be com¬ 
pared with k =0.944. For the O’*" resonance state, the momentum A;=0.086; 
and the density distribution peaks at 77 ~ 1.1. 
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Figure 14: Phase portraits for the 2+ (left) and 4+ (right) resonance states in ^Be. 


6.3. Classical regime 


In this subsection, we consider the high energy excited state of the two- 
cluster systems. So far we considered excited states with the energy less 
than 10 MeV. Now we will look for the range of energy, where phase portraits 
have an evident ’’classical shape” or where all quantum trajectories approach 
classical trajectories. 

Let us consider the excited O’*" states of ^Be with energy £^=60.21 MeV, 
Figure [T^ This energy corresponds to the dimensionless momentum k = 


^^=2.34. One can see that the maximum of the density distribution is 
observed for such value of rj which is very close to the corresponding value of 
momentum k. This means that the classical regime is realized for a relatively 
small value of the excitation energy of the two-cluster system. Besides, one 
can see that the density distribution, presented in Figure 15, is similar to 


that for free motion of a particle with a large value of the momentum k (cf. 
with Fig. and eq. (pT|). It is important to underline, that the classical 
regime becomes valid for a moderate value of the coordinate 

The correspondence between a quantum density distribution and its clas¬ 
sical limit is observed for other values of the total angular momenta in the 
^Be and for all angular momenta in other nuclei as well. The larger is the en¬ 
ergy of the excited state, the closer is the density distribution to the classical 
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trajectory. 
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Figure 15: Phase portrait for O’*' excited state in ®i?e with energy i?=60.21 MeV and 
dimensionless momentum fc=2.34 (this value is indicated by arrow). 


7. Conclusions 

We have studied the trajectories of two-cluster nuclear systems in phase 
space. The Bargmann-Segal transformation has been used to map wave func¬ 
tions of two-cluster systems in the coordinate space into the Fock-Bargmann 
space. The two-cluster systems have been studied within a microscopic 



0.000 

0.005 

0.010 

0.015 

0.020 

0.025 


32 























model which makes use of a full set of oscillator functions to expand a wave 
function of relative motion of interacting clusters. The dominant two-cluster 
partition of each nuclei was taken into consideration. The input parame¬ 
ters of the model and nucleon-nucleon potential were selected to optimize 
description of the internal structure of clusters and to reproduce position of 
the ground state with respect to the two-cluster threshold. We have consid¬ 
ered a wide range of excitation energies of compound systems, but special 
attention was devoted to the bound and resonance states. It was shown that 
bound states and narrow resonance states realize themselves in a very com¬ 
pact area of the phase space. Phase portraits of the excited states with large 
value of the energy have a maximal value along the line which coincides with 
a classical trajectory. 

We have also considered two model problems, a harmonic oscillator and 
free motion in a three-dimensional case, which helped us to understand the 
dynamics of real physical systems in the phase space. 
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